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Convergence Rates 
in Homogenization of Stokes Systems 

Shu Gu* 


Abstract 

This paper studies the convergence rates in and of Dirichelt problems for 
Stokes systems with rapidly oscillating periodic coefficients, without any regularity 
assumptions on the coefficients. 
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1 Introduction and Main Results 


The purpose of this paper is to study the convergence rates of Dirichlet problems for 
Stokes systems with rapidly oscillating periodic coefficients. More precisely, we consider 
the following Dirichlet problem for Stokes systems associated with matrix A, 


1 

( G(w) + Vpe = F 

in D, 



div Ue = g 

in D, 

(1.1) 

1 

[ u, = f 

on dQ, 


with the compatibility condition 



[ 9- [ f-n 

= 0, 

(1.2) 

Jn Jan 




where n denotes the outward unit normal to dVt and D C is a bounded domain. We 


note that the Dirichlet problem (11.11) is used in the modeling of flows in porous media. 
Here £ > 0 is a small parameter and the operator is dehned by 


G = —div(H(a;/£)V) 
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A 

dxi 



d 

dxj 


(1.3) 
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with 1 < i, j, a,/9 < (i (the summation convention is used throughout). We will assume 
that the coefficient matrix A{y) = {<\j {y)) is real, bounded measurable, and satishes the 
ellipticity condition: 

< -kt for !/ 6 R" and { = «”) 6 (1.4) 

/i 

where y > 0. We also assume that A{y) satishes the periodicity condition, 

A{y + z) = A{y) for ?/ G and z G (1.5) 

No symmetry condition on A{y) is needed. A function satisfying fll.5l) will be called 
1-periodic. 

By the homogenization theory of Stokes systems (see [21E]), under suitable conditions 
on F, f and g, it is known that 

Us Uq weakly in R'^) and Pe — j Ps Po — j Po weakly in h^(f2), 

Jn Jn 

where (mo,Po) ^ R'^) x L^(f2) is the weak solution of the homogenized problem with 

constant coefficients, 

( Co{uo) + Vpo = F in fl, 

< div Uq = g in f2, (1.6) 

\ Uo = f on Oil. 

The primary purpose of this paper is to investigate the rate of convergence of Hue—^01^2(0), 
as e —)■ 0. The following is the main result of the paper. 

Theorem 1.1. Let Ll be a bounded domain. Suppose that A satisfies the ellipticity 
condition ( [i.^D and periodicity condition Given g G H^{Q) and f G 

satisfying the compatibility condition U.S\) . for F G L^(f2;R'^), let {ue,Ps), (uo,Po) be weak 
solutions of Dirichlet problems M.Oi) . respectively. Then 

llw ~ '^^o||l2(o) < C'£||'tto||ir2(Q), (1.7) 

where the constant C depends only on d, p, and 12. 

Theorem 11.11 gives the optimal 0{e) convergence rate for the inverses of the Stokes 
operators in operator norm. Indeed, let '. F E T^(12) Ue, where T^(f2) = {F G 
L^(f2;R‘^) : div(F) = 0 in 12}, and denotes the solution of (II.ip with F G L^(12;R'^) 
and g = 0, f = 0. Then it follows from fll.7p and the estimate ||uo||h 2 (o) < C'||F|| 2 , 2 (q) 
that 

\\Ts - Zo||L 2 (Q)^i 2 (Q) < Ce, 

where Tq : F ^ Uq. 

In this paper we also obtain 0{y/e) rates for a two-scale expansion of (u^, p^) in 
Let (x, tt) denote the correctors associated with A, dehned by fl2.5l) . and S'g the Steklov 
smoothing operator dehned by fl2.ip . 
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Theorem 1.2. Let fl be a bounded domain. Suppose that A satisfies ^1-4^ and U.5\) . 
Let {us,Pe) and {uo,po) be the same as in Theorem M . 1\ Then 


Wue — Uq — £X'^*S'£(VMo)||_ffi(r2) < C'\/£||mo||//2(q), (1.8) 

where and Uq is the extension of Uq defined as in H3.1\) . Moreover, if 

InPe = InPo = 0. then 

\\Pe - Po - [T^^SfiVno) - J 7r^SfiVuo)^\\L2(Q) < Cy/e\\uo\\m{n), (1-9) 

where M{x) = n^x/e). The constants C in U.Si) and / f j..9l) depend only on d, /i, and 12. 

We now describe the known L^ convergence results on Dirichlet problems for general 
elliptic equations and systems with rapidly oscillating periodic coefficients. Consider the 
Dirichlet problem for the scalar elliptic equation Cfiufi) = —div(74(a;/£)VM£) = F in a 
Lipschitz domain 12 with = f on 512. It is well known that 

\\ue — uo\\L‘ 2 (n) < Ce {||V^Mo||L 2 (f 7 ) + II VMo||L°°(ar 2 )} • (1-10) 

To see (I1.10|h one considers the difference between Ue and its hrst order approximation 
Mo + £X^Vmo and let 

Ve = Ue-Uo-ex"'^Uo. ( 1 . 11 ) 

To correct the boundary data, one further introduces a function w^, where is the 
solution to the Dirichlet problem: Cfiwfi) = 0 in 12 and = —ex^'Vuo on 512. Using 
energy estimates, one may show that ||m£ — Ws\\h^(q) < V^mo||l2(o). The estimate 

(ll.inj) follows from this and the estimate ||ry£||L°°(o) < CellVMo||L“(ao), which is obtained 
by the maximum principle (see e.g. m)- More recently, Griso Hi was able to establish 
the much sharper estimate fll.7l) . using the method of periodic unfolding. We mention 
that in the case of scalar elliptic equations with bounded measurable coefficients, one may 
also prove fll.Tp by using the so-called Dirichlet corrector. In fact, it was shown in [9] that 

llw - Mo - - a:|vMo||iri(o) < Ce\\uo\\H2{n), (1-12) 

where $£( 0 ;) is the solution of £e(<he) = 0 in 12 with $£ = a; on 512. In the case of elliptic 
systems, the estimates (I1.12|) and thus (II.7|) continue to hold under the additional assump¬ 
tion that A is Holder continuous. Moreover, if A is Holder continuous and symmetric, it 
was proved in [8] that 

ll^e|lrfi/2(n) < Ce\\uo\\H^(n). (1-13) 

The approaches used in Pi rely on the uniform regularity estimates established 
in mm and do not apply to operators with bounded measurable coefficients. Recently, 
by using the Steklov smoothing operator, T.A. Suslina was able to establish the 

0{e) estimate (II.7p in for a boarder class of elliptic operators, which, in particular, 
contains the elliptic systems in divergence form with coefficients satisfying the elliptic- 
ity condition afffffj > /i|^P for any ^ G Since the correctors y may not 
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be bounded in the case of nonsmooth coefficients, the idea is to consider the two-scale 
expansion 

Ve = Ue-UQ-eX^Se(yUQ), (1.14) 

where Se is a smoothing operator at scale £ and uq an extension of uq to (also see 
[minus] and their references on the use of Se in homogenization). This reduces the 
problem to the control of the norm of We, where We is the solution to the Dirichlet 
problem: Ce{we) = 0 in hi and We = —ex^SeViuo) on dVL. Next, one considers 

he = We- eX^OsSei'^Uo), 

where 6e is a cutoff function supported in an e neighborhood of dQ. Note that = 0 on 
dfl and Ce{he) is supported in an e neighborhood of dfl. This allows one to approximate 
he in the norm by ho, using an 0{y/e) estimate in and a duality argument, where 
^o{ho) = jCe{he) in hi and ho = 0 on dfl. Finally, one estimates the norm of ho by 
another duality argument. 

In this paper we extend the approach of Suslina to the case of Stokes systems, which 
do not £t the standard framework of second-order elliptic systems in divergence form. As 
expected in the study of Stokes or Navies-Stokes systems, the main difficulty is caused by 
the pressure term Pe- By carefully analyzing the systems for the correctors (x,vr) as well 
as their dual (</>fc£, q^j) (see Lemmas 3.1 and 3.3), we are able to establish the 0(^/e) error 
estimates, given in Theorem 11.21 for the two-scale expansions of {ue,Pe) in x L^. This 
allows us to use the idea of boundary cutoff and duality argument in a manner similar to 
that in [T3] . 

The paper is organized as follows. In Section 2 we recall a few basic properties of 
the Steklov smoothing operator Se as well as the homogenization theory for Stokes sys¬ 
tems with periodic coefficients. In Section 3 we study uo + ex^SeVuo as the hrst order 
approximation of Ue- We introduce the dual correctors (<F, q) and use energy estimates 
to establish the estimate fll.Sp in . In Section 4 we study the convergence of Pe and 
prove the error estimate fll.9p for the two-scale expansion of the pressure term. Finally, 
our main theorem Theorem 11.11 is proved in Section 5. This is done by using the idea of 
boundary cutoff and duality, and by applying error estimates obtained in Sections 3 and 
4 to the adjoint systems. 

Throughout this paper, we denote Y = [0, l)'^ and the average of / over the set E 



We will use C to denote constants that may depend on d, fi, or hi, but never on e. 
Acknowledgement. The author would like to thank referees for their very helpful 
comments and suggestions. 
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Preliminaries 


2.1 Smoothing in Steklov’s sense 

Let S'e be the operator on L‘^(W^) given by 

{Seu){x) = u{x — ez)dz (2.1) 

and called the Steklov smoothing operator. Note that 

ll>S'£M||L2(Kd) < ||M||i2(]Kd)- 

Obviously, D'^SsU = S^D^u for u G H^(W^) and any multi-index a such that |a| < s. 
Therefore, 

The following are a few properties of Steklov’s operator; see [T^fH] . 

Proposition 2.1. For any u G we have 

IIS'eM — M||£,2(Rd) < Cell VM||i;,2(Rd), 

where C depends only on d. 

We will use the notation f^{x) = f{x/e). 

Proposition 2.2. Let f{x) be a 1-periodic function inW^ such that f G Lf{Y). Then for 
any u G 

\\f^^eU\\L^(Rd'^ < ||/||L2(y)||M||2^2(Rd). 


2.2 Homogenization of Stokes systems 


We refer the reader to [2116] for details of weak solutions and homogenization theory of 
Stokes system. 

Let be a bounded Lipschitz domain in For u,v G we dehne the 

bilinear form as(-, ■) by 


a^(u,v) 



du^ dv°‘ 
dxj dxi 


dx. 


For F E H and g G we say that (u^jPe) G x L^(f2) is a weak 

solution of the following Stokes system in f2. 


if for any ip G W^), 


Ce{Ue) + Vpe = F 
div = g, 


ae{Ue,^) - Pe div(v3) = {F, 


and div(M£) = in (in the sense of distribution). 


( 2 . 2 ) 
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Theorem 2.3. Let fl be a hounded Lipschitz domain in Suppose A{y) satisfies the 
ellipticity condition Let F G g G L‘^{Lt) and f G M'^) satisfy 

the compatibility condition Then there exist a unique G and Pe G 

L'^ifil), unique up to constants, such that {ue,Pe) is a weak solution of H2.^) and = f 
on dQ. Moreover, 


+ WP^ ~ £ P.lli>(n)<C'{||F|| H-\n) + WaWi^n) 

where C depends only on d, pi, and f2. 


+ 


(2.3) 


Theorem 12.31 is proved by using the Lax-Milgram Theorem. We mention that if 
is and A is a constant matrix, the weak solution {u,p), given by Theorem 12.31 is 
in x H\Ll), provided that F G L‘^{Ll;R^), g G and / G 

Moreover, 

ll'*^l|r/2(n) + ||Vp||L2(f^) < c|||F||i2(Q) + ||5'||//1 (q) + ||/||/f3/2(aQ)|, (2.4) 

where C depends only on d, p, and hi (see e.g. 0). 

We denote by the closure in of C:^fiY;M.‘^), the set of C°° 

1-periodic and M'^-valued functions in Let 


®per (aP y 0) 

where ; R'^). Define 


ay, 

o.. o.. dy, 


>Y 


dxj dxi 


Dper(T) = ju G Hl^fiY-Rh^) : div(M) = 0 in F and ^ n = o|. 

By applying the Lax-Milgram Theorem to aper(^, fi) on the Hilbert space Lper(F), one may 
show that for each 1 < j, /? < d, there exist 1-periodic functions {Xj, vr^) G ^ 

L^jj^(R'^), which are called the correctors for the Stokes system fl2.2p . such that 

^ ) + Vtt^ = 0 

< div = 0 

1^71^ = 0, J^Xj =0, 

where {y) = = yfit), - ■ ■ , 1, • • • ,0) with 1 in the position. Note that 

\\x^j\\m{Y) + lkf||L2(y) < C, 

where C depends only on d and /i. The homogenized system for the Stokes system (12.21) 
is given by 

'^o(mo) + Vpo = -^ , . 

div uo = g, 


m 


m 


(2.5) 
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where Cq = —div(74V) is a second-order elliptic operator with constant coefficients, and 
= (^j), with 

sS'' = “p-(h + d.x? + ^’“). 


We remark that (A)* = A*, and the effective matrix A satisfies the ellipticity condition 

for any ^ G and /xi depends only on d and /i. The following 
is a homogenization theorem for the Stokes system. 

Theorem 2.4. Suppose that A{y) satisfies ellipticity condition o and periodicity con¬ 
dition M.^) . Let Q be a bounded Lipschitz domain. Let {us,Ps) G x L^(f2) be a 

weak solution of Hl.l\) . where F e g e L^(f2) and f G Assume 

that J^Pe = 0. Then, as e ^ 0, 

Ue Uq strongly in L‘^{Q;Mfi), 

Ue Uq weakly in Mfi), 

Pe po weakly in 

A{x/e)'Vue AVuq weakly in 

Moreover, {uq,Pq) is the weak solution of the homogenized problem 


3 Convergence rates for Ue in H 


1 


From now on we will assume that is a bounded domain with boundary of class 
F G g G H^{Lt), and / G We fix a linear continuous extension 

operator 


and let 

so that Mo = Mo in ^ and 


Wo — EqUq, 


l|wo||M2(Rd) < C'||wo||iy2(n), 

where C depends on f2. We introduce a first order approximation of m^, 

Ve = uo + ex^SfiWuo). 

Let G x be a weak solution of 

Cfiwe) + Vr^ = 0 in hi, 

div(t(;£) = £ div(x*^5£VMo) in hi, 

We = £X^S'e(VMo) On dkl. 


(3.1) 

(3.2) 
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(3.3) 







We will use Wg to approximate the difference between and its first order approximation 
Ve- To this end, for 1 < i,j, a, /3 < d, we let 


d 


itiv) = <(!/) + Uf) - SS'’ 


(3.4) 


Note that 6“^ is 1-periodic. By the definition of y and A, 6"^ G L‘^{Y) satisfies 


Kfiy) dy = 0 - 


and, for each 1 < q;,/5,j < d, 


A 

dVi 





A 

^ dVi 

_ 

dyi 


A(l/) 


dxf 

dyk 



dPf\ 
dyk ) 


+ 



Lemma 3.1. There exist G Hp^^{Y) and G Hp^^{Y) such that 


b-f = A(4g) + 3 (,s 




dyk 


dya 


ijj 


a.nd <. = -4g. 


(3.5) 


(3.6) 


Moreover, 

W^kijWL^iY) + lkt§IU2(y) < C, (3.7) 

where C depends only on d and pi. 

Proof. Fix 1 < i,j,/3 < d. There exist = (A) G H^^^(Y-,R^) and G H^ev(X) 
satisfying the following Stokes system, 

^ A/^ + Vgg = fej in F, 

< div(/^) = 0 in F, 

[ fijdy = 0, 

.Jy 


where — ibY). We now define 


(3.8) 



Clearly, 6 HUY) and = -4“/,. Note that, by ga and 631, ^ (F) 


satisfies 


■ kij 


per \ 


■ kij 


■ ikj 


dVi 


per \ 



d (dq 


dya V dyi 


+ 


8Kf 


d 


dq 


dyi dye 


TTl- 


dyi 


(3.9) 


dya y dyi 

It follows by the energy estimates that 


dfTf 




is constant. Hence, 


® (A?)) = Kf - ^( 4 )- 


dyk^^^^^’ dykdyk^''^ ’ dy^dyk^'^^ , 

Furthermore, since ||Xj ||Hi(y) ^ C, then 

\\^Hj\\L^{Y) + WQijWh'^iY) < C\\lq^\\u(j) < C, 
where C depends only on d and y. This completes the proof. 


□ 


B B ^ 

Remark 3.2. Recall that and are both 1-periodic. By and the fact that 


Bq' 

is constant, we see that tt^ and differ only by a constant. Since fy tij = 0, we obtain 
the following relation, 


g ^4 

TT • = - - 

’ Sy, 


(3.10) 


Lemma 3.3. Let Q be a bounded domain. Suppose that A satisfies ellipticity con¬ 
dition a and periodicity condition / (i.5]) . Given g G and f G M'^) 

satisfying the compatibility condition / li.dl) . for F G let {ue,Pe), {uq,Pq) and 

{ws,Te) be weak solutions of Dirichlet problems U.l\) . U.6\) and H3. 3)) . respectively. Then, 


\\Ue - Mo - < C'^IImqII ^2(0), 

where C depends only on d, y, and H. 

Proof. Let 


(3.11) 


= u. 


Then 


uo- ex^SfiVuo) Fwe- 
div(.^£) =0 in H and = 0 on dfl. 
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Now we compute 


{C-e{z,)T = 


d 


d[Pe - PO + n] _ 

dxi 


dXa 

d r 


'^ajS olB ( I \ 

“a - %(xle) 


i^g 

dXn 




dxk 

d[Pe -Po + Te] _ 

dxi 


du^ 


dxi 






dXr 


OlB / / \ 

ap - «p 


du^o _ g du^o 


dxi 


dxn 


d 


du^c 


d 


+ ^ bf(x/e)S,^ +.— al-ixMxfixMS, 


dxi \ dxj J dxi 

Using Lemma [3.11 we may write 


O'^Uq 

dxkdxj 


'dxkdxj 


d 


du^ 


d 


^ I I = -7^ 


dxi 


dx 

= h + h. 


A. 

dxi 




S', 


du^ 

dxj I (3.12) 


Since that 


ikj ’ 


ii = 


02 




dxidxk 

d 


dup^ 

dXn 


■A 

dxi 




= —e 


dxi 




'dxjdxk 


For the second term in the RHS of fl3.12p . we have 


h = 


d d 


dxn \ dxi 




dx^^ 

= 4 - A I ^ 


— ieQHxle)S ^'^0 

dxi 1 ^ dxadxj 


A 

dxi 


dXadXj 


(3.13) 


In view of fl3.10l) . for the hrst term on the RHS of fl3.13l) . we obtain 


h = 


A 

dXr 




+ 


A 

dXr 


/3, , ac. 

eq^Ax/e)S, 


dxjdxi 


(3.14) 
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Putting altogether, we have shown that 




A 

dXr, 




Pe-Po- Ti^Ax/e)Se'^ - eqAx/e)Se ^ A + n 




dxi 


' dxjdxi 


= e 


A 

dxi 


a^^{xle)xT{xle) - ^%{xle) 
d‘^uf 


S 


'dxjdxk 


-» 


(3.15) 


dXadXj 


A 

dxi 


A “ A (^/^) 


du^ c ^^0 


Since 2 ;e G and div( 2 ;£) = 0 in fi, it follows from fld.lSp by the energy estimate 

fl2.3p that 


c / \Vzs\‘^dx < / [|x(a:/e:)| + |<h(x/e:)|]S'£(V^Mo) 


dx 




q{x/e)Se(y‘^Uo) 


dx + 


Vmo - ^^(Vmo) 


dx. 


Now we apply Propositions 12.1112.21 as well as fl3.2p . This gives 

II V 2 :e||i, 2 (Q) < Ce (||x||L 2 (y) + ||<h||i 2 (y) + ||g||i 2 (y) + l) || V^Mq || L 2 (Rd) 

— e'en V^mo||l2(r'^) 

< Ce||Mo||H2(o), 

where C depends only on d, p and 12. Hence we have proved the desired result, ||2:e||Hi(o) 
Ce||Mo||H 2 (Q), and completed the proof. 

For r > 0, let 

= {x G : dist(x, 912) < r}, 

12j, = {x G 12 : dist(x, 912) < r}. 

We choose two cut-off functions 9^{x) and 9^{x) in satisfying the following conditions. 


0^ G C'“(R‘^), supp(6'£) c (912)^, 0 < e^{x) < 1, 

9e\dQ. = 1, iVdel < n/e, 


(3.16) 


and 


6e G C“(R'^), snpp(6'£) C (912)2e, 0 < ^^(x) < 1, 

^^(x) = 1 for X G (912)£, iVdgj < 'k/e. 

The following is an estimate for integrals near the boundary, see [14] for example 


(3.17) 
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Lemma 3.4. Let fl G be a hounded domain. Then, for any function u G H^{VL), 


/ \u\^dx <Cr\\u\\H^(a)\\u\\L^(a)■ 

Moreover, for any 1-periodic function f G L^{Y) and u G 


\f^\ IS'gWl dx < C'e||/||i 2 (y)||M||j^i(Rd)||M||^ 2 (Rd), 


J{dU)2e 

where C depends only on f2. 

We are now ready to give the proof of fll.Sp . 


Proof of estimate 1^1.8\) . By Lemma [3.31 the problem has been reduced to estimating 
Ws in Notice that by the energy estimate (I2.3p . 

Ikellnpo) < Ce\\x^SeVuo\\HG2(9n) + C£||div(x^S'£VMo)||L2(n) 

< C^II^'eX^S'eVrtollHpO) + C'£|lX^VS'e(VMo)||L2(f^) 

< Ce{\\x^S,{Vuo)\\LHn) + ll(V0,)x"^e(VMo)||L2(o) 

+ e-^\\9,{VxrS,{Vuo)\\LHn) + \\x^S,{V^uo)\\mn) 

< Ce|||Mo||_H-2(Rd) +e ^||x^5'£(Vmo)||l2(q^) +£ ^||(Vx)'^5'£(Vmo)||l2(Oj) 


(3.18) 


where we have used Proposition 12.21 for the fourth inequality and Lemma [3.41 for the last. 
We point out that the fact div(x) = 0 in is also used for the second inequality in 
03.181) . Therefore, 

\\Ue — Uq — eX^Se (Vmo) 11^1(0) < Ikellr/pO) + 

< C\fe\\uo\\H2{n), 


where C depends only on d, /i, and fl. This completes the proof. 


□ 


4 Convergence rates for the pressure term 

To prove estimate Ol.bh . we hrst recall that if {ue,Pe) £ x is a weak 

solution of the Stokes system 01.Ih . then 


(4.1) 


WPs — T PeWl^in) < C||Vpe||_H--i(Q) < C< ||F||_H--i(0) + ||m. 




where C depends only on d, /i, and fl (see e.g. USD- 
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Proof of estimate Since = JqPo = 0, using fl4.ip and fl3.15l) . we see that 


\\Pe -Po- {TT^SeVuo + eq^SeV'^Uo -Te) - -I {iT^SeVuo + Sq^SeV'^Uo - T^) 


|L2(Q) 


r. 




< C|| V Pe-Po- (Vmo) - eq^Se (V^Mo) + 

< C'{||Vz,|U2(n) +HI(lx1 + m + (V'«o)|L 2 (f,) + \\Ss{Vuo) - VmoIIl^cq) 

< Ce\\uo\\H2{n), 

(4.2) 

where the last inequality follows from the proof of Lemma [3.31 Note that by Propostion 

Oand (I32D, 


ellg^'S'gV Mo - f Mo||L2(n) < C£||Mo||ir2(Rd) < Cs\\uo\\H^n). (4.3) 

Jq 

Also, by the dehnition of {w^,Te) and fl4.1l) . 


:||l2(o) < C'll Vrell^-po) < C\\'\/We\\L^(n) < C\/e\\uo\\H2{n), (4.4) 


r. - r. 


where the last inequality follows from (I3.18p . By combining (14. 2 h , (14. 3 h and (14.4p , we have 
proved that 


TT 


WPs -Po- 
This completes the proof. 


‘S'e (VMq) “ T (VMq) ||L2(n) < C'V^||mo||_H'2(q). 

Jn 


□ 


5 Convergence rates for in 

To establish the sharp 0{e) rate for in in view of (I3.11h . we obtain 
ll^e — Mo — ^X^SeVuo + We\\L^(n) < Ce\\uo\\ h'^( n)■ 


Using Proposition 12.21 and fl3.2p . 

||x'^S'£VMo||L2(f^) < C||x||l,2(y)|| VMo||L2(Rd) < C'||mo||// 2 {n) • 


Thus, 


Me — Mo||L2(f2) < Ce\\uo\\H2(n) + ||m;£||i,2(q). 


and it remains to estimate ||M;£||i 2 (Q). 


(5.1) 


Lemma 5.1. Let Q be a bounded domain. Suppose that A satisfies ellipticity con¬ 
dition fjl.fl) and periodicity condition M.^) . Given g G H^iVt) and f G 
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satisfying the compatibility condition M.2^) . for F G let {u^,Pe), {uo,Po) be 

weak solutions of the Dirichlet problems M.b\) . respectively. Then 


\Ue — Uo — £(1 — 6e)x^Se (Vmq) \\h^{Q) < C ^/e\\uo\\ 


(5.2) 


and 

WPs -Po- (1 - 0,)7r^Ss (Vmo) - / (Vmo) ||l2(q) < C^/e\\uo\\H2(n), (5.3) 

where C depends only on d, fv, and f2. 

Proof. Note that 

\\£d^X^Se'Vuo\\m(n) < C£||x^*S'£(Vmo)||l2(q2j) + C£||x^5'e(V^Mo)||L2(Q2e) 

+ C^II(Ix^I + I(Vx)^I)^s(V«o)||l^(o..) (5.4) 

< Cy/e\\uo\\H2{n), 


where we have used Lemma [33] and Proposition fl2.2l) for the last inequality. This, together 
with estimate fll.Sp . gives fl5.2p . 

Similarly, using Lemma 13.41 we see that 

\\e,7r^S,Vuo\\l^n) <C [ |7r^S'3VMo)r < C'^I|wo|Ih2(o)- 

d{dn) 2 e 


This, together with estimate (II.9p . gives 
Proof of Theorem \1.1[ In view of fl5.ip , it suffices to show that 

||'?i'e||L2(Q) < Ce\\uo\\H^{n)- 

Furthermore, let 

(pe = e9eX^SeVuo- 

Since ||0e||i,2(Q) < C'e||Mo||H 2 (o), it is enough to show that 

||h£||L2(o) < Ce\\uo\\H2(n), 


□ 


(5.5) 


where — (ps- 

To this end, we hrst note that by the dehnition of {ws,t^) in fl3.3p . the functions 
ir]£,T^) G iLo(f2;M'^) x L^(ff) satisfy 


k^eiVe) + Vr^ = -C^(pe m n, 

div pe = £ div ((1 - 6s)x^SsVuo) in kl, 
rjs = ^ oil dkl. 


(5.6) 
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Let (r^o, To) G Hq{Q] M.^) x be a weak solution of the homogenized Dirichlet problem 


( ^oiVo) + Vto = in n, 

< div r]o = e div ((1 - 9s)x^SsVuo) in hi, 

[ rjo = 0 on d^l 


(5.7) 


To estimate rj^ — rjo, we consider the following duality problems. For any H G L^(hl;R‘^), 
let {pe,ae) G iLo(hl;M'^) x L^(hl) be the weak solution of 


(Cl{p,) + Vae = H infi, 

< div ps = 0 in hi, 

[ Pe = 0 on c?hl, 

and (pojCTo) G (i/^(hl;M'^) fl iLo(hl;M^)) x iL^(hl) the weak solution of 

f £o(Po) + Vuo = H in fi, 

< div po = 0 in hi, 

I Po = 0 on dQ, 


(5.8) 


(5.9) 


with 




CTo = 0. 


Here we have used the notation; £* = —div(H*(a;/e)V) and £q = —div(y4*V). We note 
that Lemma [5.11 continues to hold for £*, as A* satisfies the same conditions as A. Also, 
by the estimates fl2.4l) for Stokes systems with constant coefficients in domains. 


||po||ir2(n) + ||(To||r/i(o) < C 

As a result, we have 

||Pe ~ Po ~ ^(1 ~ (Vpo) 11^1(0) < C'\/i||po||H 2 ( 0 ) < C y/e\\ H \\^2 , 

and 




where (x*,7r*) denotes the correctors associated with the adjoint matrix A* 
Let T = and 

r = div (£(1 - e^)x‘'S^Vuo)- 
Note that by fl5.6l) . (15.71) . (15.81) and (15.9p . 


(5.10) 


- ao - [(1 - (V^) - / (Vpo) 11| L2(0) ^ Cy/e\\H\\L2(^Q'j, (5.11) 

L Jo, 


/ 77 • (Pe — Po) — (^5 Pe ~ Po)_ff-l( 0 ;Rd)xml(f^;®‘*) “ / r(<T£ — (To) /r 
'n Jn (5-12j 

= A + A- 
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For the first term of the RHS of fl5.12p . because e H ^(r2;R‘^) is supported in 
and 1 — 6'e = 0 in we obtain 

Jl = Pe — PO — ^(1 ~ Po)) 


Therefore, 

I'^ll < ll^llir-i(O) lIPe ~ Po ~ ^(1 “ ^e)X*^Se (Vpo) IliTpO) 

< C\\e9eX^Se'Vuo\\Hi{n) ||l2(q) (5.13) 

< Ce\\uo\\H2{n) ||iir||L2(Q) 

where the second inequality follows from fl5.10l) . and the last inequality follows from the 
analog of (15.dh (with 9s replaced by 9s). For the second term of the RHS of (I5.12p . we 
recall that div (y) = 0. Hence, 


" - + e(i - . r.. r„ 

Since J^T = 0, for any constant E, 


J 2 = - / F(cr£ -ao + E) = - / [Fi + F2](cr£ - cxo + E). 

Jn Jn 


We split J 2 as two integrals, for the first integral, again since 1 — Pg = 0 in (c?n)e and Fi 
is supported in (c?H)£, just as we did for Ji, 



CTo + E) 



{as-ao-{l-9s)T^*^Ss{yPo) + E). 


Now, if we choose the constant E as E = j^'n*^Ss (Vpo), then 


Fi((j£ — aQ + E) 


Fi a, 


- no - [(1 - 9s)7^*^Ss (Vpo) - / (Vpo) ’ 

L Jn ^ 

< C||Fi||L2((ao),) \fe\\H\\L2{n) 

< C(\/i||y||L2(Y)||V'Uo||Hi(K'i)) (\/i||Ff ||l2(o)) 

< Ce\\uQ\\H‘2(n)\\H\\L‘^(n)i 

where we have used (15.lip and Lemma [3.41 For the second integral in J 2 , we have 


(5.14) 


/ -<To + E) 

Jn 

< l|r2||L2(0)||c’'£ — <70 + J '^*^Ss (VPo) ||l2(0) 

Jn 


(5.15) 


< Ce\\uo\\H2{n)\\H\\L2(n), 
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where for the last inequality we have used 


\\<^e ~ <^0 + T (Vpo) ||L 2 (n) < HcTg ||+ ||cro||L 2 (Q) + II T (Vpo) ||L 2 (n) 
Jo. Jo. 

Therefore, by combining fl5.13p - fl5.15p . we have proved 


H{Ve - Vo) 


< Ce\\uo\\HVn)\\H\\L2(n) for any H E 


By duality this implies that 


lihe ~Vo\\L2{n) < Ce\\uo\\H2(n)- 


(5.16) 


(6.17) 


Finally, the problem has been reduced to the estimate of ||?7 o||l2(o)- This will be done 
by another duality argument. Let (po,o'o) be dehned by fl5.9p . Then 


[ H-po 

= 

Jn 



Po)ir-i(0;R‘i)xr/i(0;R‘') “ / Tao 
LiCTo 


^ K^)Po)_H'-l(0;Rd)xHl(0;R‘*)l + 

= K, + K2 + iLs, 


^2(^0 


(5.18) 


where \l',r,ri and T^are as denoted above. Notice that again by Lemma 13.41 and the 
analog of fl5.4p (with 9^ replaced by 9s), we have 


Ki < ||4/||iy-i(n)||po||r/i(Oe) 

< C\\s9sX^SsVuo\\Hi{n)'/£\\po\\H2{Q) 

< C(\/i||no||r/2(n))(\/i||po||rf2(o)) 

< Ce||Mo||ir2(Q)||i/||2,2(n). 


(5.19) 


Similarly, again by Lemma El 


^2 < ||ri||L2((aQ)^)||ao||L2(n^) 

< C(\/£|| x||L 2 (y) 11^0II H 2 (Kd))(\/i||cro 11^1(0))) 

< C£||Mo||iy2(n)||i/||L2(f^), 


and 


K 3 < ||r2||L2(Q)||cro||L2(f^) < C£||Mo||ir2(n)||iL||i2(Q). 

By combining fl5.19l) - fl5.21l) . we obtain 


(5.20) 


(5.21) 


H-Vo 


< C£||Mo||iy2(Q)||iL||L2(Q), 
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which, by duality, leads to 


||^o||l2(q) < Ce||Mo||//2(o)- 


(5.22) 


Hence we have proved that 

||'M^£||l2(q) < ||?7e — '/7o||L2(n) + ||^7o||L2(n) + ||0 £||l2(o) < Ce\\uQ\\H'2{yi). (5.23) 

The proof is hnished. □ 
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